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Abstract 

Substantial collective flow is observed in collisions between large 
nuclei at high energy, and the data are well-reproduced by perfect fluid 
dynamics. In a separate development, calculation of the dimensionless 
ratio of shear viscosity rj to entropy density within AdS/CFT yields 
ri/s = l/47r, which has been conjectured to be a lower bound for any 
physical system. It is shown that the transition from hadrons to quarks 
and gluons has behavior similar to helium, nitrogen, and water at and 
near their phase transitions in the ratio r]/s. Conversely, there are 
indications that the ratio of bulk viscosity C to entropy density may 
have a maximum in the vicinity of the phase transition. Therefore it 
is possible that experimental measurements can pinpoint the location 
of the transition or rapid crossover in QCD via the ratios rj/s and C/s 
in addition to and independently of the equation of state. 
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1 Introduction 



A transition from a phase of hadrons to a phase of quarks and gluons with 
increasing temperature and/or baryon chemical potential has been studied 
theoretically for three decades [1]. Calculations with effective hadronic field 
theories and with perturbative QCD give consistent and intuitively reason- 
able estimates of where this transition occurs, but at the present time only 
lattice gauge theory calculations on lattices of finite size can yield quan- 
titatively accurate numbers. For two fiavors of massless quarks the phase 
transition is second order. For three fiavors of massless quarks the phase 
transition is first order. For the real world with nonzero masses for up, down 
and strange quarks the answer is not definitively known yet, but the answer 
is likely to be a rapid crossover from one phase to another without a rigorous 
thermodynamic phase transition, at least for zero baryon chemical potential 
[2]. Indeed, there may be a line of first order phase transition of the plane of 
temperature T versus baryon chemical potential /j, starting from the chemi- 
cal potential axis and terminating at some critical point in the T — /i plane 
[3, 4, 5]. The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory was constructed explicitly to create quark gluon plasma. After 
nine years of operation, what have the experiments told us? 

One of the amazing discoveries of experimental measurements of gold on 
gold collisions at RHIC is the surprising amount of collective flow exhibited 
by the outgoing hadrons. Collective flow is evidenced in both the single- 
particle transverse momentum distribution [6] , commonly referred to as radial 
flow, and in the asymmetric azimuthal distribution around the beam axis 
[7], commonly referred to as elliptic fiow and characterized by a Fourier 
coefficient called V2. Elliptic flow was expected to be much smaller at RHIC 
than it was at the lower energies of the SPS (Super Proton Synchrotron) at 
CERN [8]; in fact, it is about twice as large. Various theoretical calculations 
[9] support the notion that collective fiow is mostly generated early in the 
nucleus-nucleus collision and is present before partons coalesce or fragment 
into hadrons. Theoretical calculations including only two-body interactions 
between partons cannot generate sufficient flow to explain the observations 
unless partonic cross sections are artificially enhanced by more than an order 
of magnitude over perturbative QCD predictions [10]. Thus quark-gluon 
matter created in these collisions is strongly interacting, unlike the type 
of weakly interacting quark-gluon plasma expected to occur at very high 
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temperatures on the basis of asymptotic freedom [11]. On the other hand, 
lattice QCD calculations yield an equation of state that differs from an ideal 
gas only by about 10% once the temperature exceeds 1.5Tc, where Tc ^ 
175 MeV is the critical or crossover temperature from quarks and gluons 
to hadrons [2]. Furthermore, perfect fluid dynamics with zero shear and 
bulk viscosities reproduces the measurements of radial flow and f 2 quite well 
up to transverse momenta of order 1.5 GeV/c [12]. Parametric fits to the 
transverse momentum spectra of hadrons, such as pions, kaons, and protons, 
result in average transverse fluid flow velocities of more than 50% the speed 
of light and chemical freezeout temperatures on the order of 160 MeV. These 
results have been interpreted as strong indicators of early thermalization and 
collective flow on a time scale of several fm/c. 

An amazing theoretical discovery was made by Kovtun, Son and Starinets 
[13]. They showed that certain special field theories, special in the sense that 
they are dual to black branes in higher space-time dimensions, have the ratio 
77/s = l/An (in units with h = kB = c = 1) where 1] is the shear viscosity 
and s is the entropy density. The connection between transport coefficients 
and gravity is intuitively clear since both involve (commutators of) the stress- 
energy-momentum tensor T^^. They conjectured that all substances have this 
value as a lower limit, and gave as examples helium, nitrogen, and water at 
pressures of 0.1 MPa, 10 MPa, and 100 MPa, respectively. Is the RHIC data, 
represented by radial and elliptic flow, telling us that the created matter has 
a very small viscosity, that it is a perfect fluid! 

The relatively good agreement between perfect fluid calculations and ex- 
perimental data for hadrons of low to medium transverse momentum at RHIC 
suggests that the viscosity is small; however, it cannot be zero. Indeed, the 
calculations within AdS/CFT suggests that rj > s/An. It will now be argued 
that sufficiently precise calculations and measurements should allow for a de- 
termination of the ratio r)/s as well as the ratio of bulk viscosity to entropy 
density (/s as functions of temperature, and that these ratios can pinpoint 
the location of the phase transition or rapid crossover from hadronic to quark 
and gluon matter. This is a different method than trying to infer the equa- 
tion of state of QCD in the form of pressure P as a function of temperature 
T or energy density e. 

The shear and bulk viscosities are rigorously given by the Kubo formulas 
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which are based on hnear response theory [14]. 

V = ^lim^/cfW'^*([7:4eeless(^),^4celess(0 (1) 

C = ^hm - fd'xe^'{[Pit,K),V{O,O)])0it) (2) 

Z w— >u (jj J 

Here T^^^ccicss represents the traceless part of the spatial components of the 
energy-momentum tensor and V = represents the trace of the tensor 
(the pressure in equihbrium). There are some interesting hmiting cases. 

Why is the entropy per baryon of the universe as large as 10^? It had orig- 
inally been presumed that the expansion of the universe was adiabatic. Al- 
though shear viscosity and thermal conductivity play no role in a Robertson- 
Walker model, bulk viscosity could. In order to address this problem, as well 
as the problem of damping of protogalactic fluctuations immediately prior to 
the period of recombination of hydrogen, Weinberg was led to study the bulk 
viscosity due to radiation [15]. Much earlier Thomas [16] had calculated the 
correction to the energy-momentum tensor due to radiation, meaning pho- 
tons, neutrinos, or gravitons, interacting with material particles with much 
smaller mean free times. Prom that result Weinberg inferred the shear and 
bulk viscosities. (He also inferred the heat conductivity which is of less in- 
terest for RHIC or LHC, the Large Hadron Collider at CERN.) For photons 

4 

C = AaT\{l-{dP/de)ny (4) 

where a = 7r^/15, Tj is the photon mean free time, and the derivative 
{dP/de)n is taken at constant baryon density. For neutrinos these results 
are multiplied by a factor of 7/8. This leads to the oft-quoted result ( — 
15(| — VgYrj. If the material particles themselves are nearly massless then 
C/rj — > 0. One must always keep in mind the conditions for which these 
results apply, namely, massless particles with a very large mean free time 
compared to material particles. In heavy ion collisions the photons escape 
from the system without interacting, so this formula is not directly relevant 
for RHIC and LHC. 

Gavin [17] derived relativistic formulas for a single species of hadrons in 
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the relaxation time approximation. 



15Ti (2.)=B^^"(^/^) 

The factor — v^)^ appears once again, but here Vg refers to the speed of 
sound of the single species of particles themselves. There is no additional ra- 
diation that provides the transport. The bulk viscosity still has the property 
that it goes to zero as — >^ |. 

Apart from entering the relativistic fluid equations, the viscosities play 
another role for first order phase transitions. They are prefactors in the 
nucleation rate [18] 

T = 1 (^\ ^(Ch + 4?7H/3)i^. _AF/r /7^ 

where AF = AnaRl/S is the free energy and i?* is the radius of a critical 
size fluctuation as given by the Laplace formula. Also a is the surface free 
energy (surface tension) and Aw is the difference in enthalpy densities of the 
two phases. This is the probability per unit volume per unit time to nucleate 
an L (low temperature) phase bubble out of the H (high temperature) phase. 
If one considers nucleating an H phase droplet in the L phase instead, one 
just needs to evaluate the correlation length and the viscosities in the L phase 
rather than the H phase. At the critical temperature, oo, and the rate 

vanishes because of the exponential. Any system must supercool at least a 
minute amount in order that the rate attains a finite value. 

It is important not only to ask what has RHIC told us about the equation 
of state, or what LHC will, but to make connections with other areas of 
physics such as cosmology and condensed matter physics. To this end section 
2 will summarize some interesting and important facts about atomic and 
molecular systems. Relativistic dissipative fluid dynamics will be reviewed 
in section 3. The shear viscosity will be studied in section 4 and the bulk 
viscosity in section 5. The relevance and advances in applying gauge/gravity 
duality methods to transport coefficients in QCD are discussed in section 
6. An overview of possible experimental observables is given in section 7. 
Concluding remarks are in section 8. 
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2 Atomic and Molecular Systems 



Shear viscosity rj is relevant for a change in shape at constant volume. Bulk 
viscosity ( is relevant for a change in volume at constant shape. Here we 
review what is known about these in the vicinity of a phase transition in 
atomic and molecular systems. 

In figures 1 through 3 we [19] plot the ratio rj/s versus temperature at 
three fixed pressures, one of them being the critical pressure (meaning that 
the curve passes through the critical point) and the other ones being larger 
and smaller, for helium, nitrogen and water. The ratio was constructed 
with data obtained from the National Institute of Standards and Technol- 
ogy (NIST) [20]. (Care must be taken to absolutely normahze the entropy 
to zero at zero temperature; we did that using data from CODATA [21].) 
The important observation [22] is that rj/s has a minimum at the critical 
point where there is a cusp. At pressures below the critical pressure there 
is a discontinuity in rj/s, and at pressures above it there is a broad smooth 
minimum. It is noteworthy that this behavior is exhibited by systems with 
very different phase transitions. In the vicinity of the critical point, helium 
is a quantum system whereas nitrogen and water are classical systems. Fur- 
thermore, nitrogen is a diatomic molecule whereas water is a polar molecule. 
It is interesting to ask whether there is a universal scaling behavior. A plot 
of {rj/ s)/{r)/ s)c versus T/Tc, shown in figure 4, indicates that there is not. 

The simplest way to understand the general behavior was presented by 
Enskog, as explained in [23]. Shear viscosity represents the ability to trans- 
port momentum. In the classical transport theory of gases 

n Tv 

~ TZfree^ (8) 

s na 

where Zfree is the mean free path and v is the mean speed. For a dilute gas the 
mean free path is large, /free ~ l/no", with n the particle number density and 
a the cross section. Hence it is easy for a particle to carry momentum over 
great distances, leading to a large shear viscosity. (This is the usual para- 
dox, that a nearly ideal classical gas has a divergent shear viscosity.) In a 
liquid there are strong correlations between neighboring atoms or molecules. 
A liquid is homogeneous on a mesoscopic scale, but on a microscopic scale 
it is a mixture of clusters and voids. The action of pushing on one atom is 
translated to the next one and so on until a whole row of atoms moves to 
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fill a void, thereby transporting momentum over a relatively large distance 
and producing a large shear viscosity. Reducing the temperature at fixed 
pressure reduces the density of voids, thereby increasing the shear viscosity. 
The shear viscosity, normalized to the entropy, is observed to be the smallest 
at or near the critical temperature, corresponding to the most difficult con- 
dition to transport momentum. This is an empirical observation. It is also 
supported by numerical calculations. For example, the shear viscosity was 
calculated for 2-dimensional Yukawa systems by Liu and Goree [24]. These 
have interesting applications to dusty plasmas and many other 2-dimensional 
condensed matter systems. Figure 5 is a plot of the shear viscosity versus the 
Coulomb coupling parameter F. Here F = Q^/aT, where Q is the charge and 
a = l/y^7m is the Wigner-Seitz radius. The shear viscosity is large at high T, 
which is a gas phase, where it is dominated by kinetic energy contributions. 
It is also large at small T, which is a Hquid phase, where it is dominated 
by potential energy contributions. The shear viscosity has a minimum near 
F = 17 where kinetic and potential contributions contribute approximately 
equal amounts. 

It has long been known that the bulk viscosity for a simple nonrela- 
tivistic gas of point particles, interacting via short-range 2-body forces, is 
vanishingly small. This result has been attributed to Maxwell [25] . However, 
the bulk viscosity is not always small compared to the shear viscosity. One 
quite relevant example is a diatomic gas. The exchange of energy between 
translational and rotational degrees of freedom during collisions results in an 
appreciable bulk viscosity, which is important in the damping of sound waves 
[23]. 

Another relevant example is a liquid-gas phase transition. The bulk vis- 
cosity can be calculated from classical molecular dynamic simulations using 
the classical limit of the Kubo formula, which reads 

C = ^ E dt J d'x{{TuiO) - (T,,)) (r,,(x, t) - (T,,))) (9) 

where Tij are the spatial components of the energy-momentum tensor, and 
the angular brackets indicate the ensemble average. This is just a measure- 
ment of the pressure fluctuations. The results of one such simulation by 
Meier, Laesecke and Kabelac [26] are shown in figure 6. The bulk viscos- 
ity (not normalized to the entropy density) shows a peak in the vicinity of. 
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though not exactly at, the critical point. The pressure fluctuations are large 
in this situation because there are two phases below the critical temperature 
with different energy densities. The molecular dynamic simulations conserve 
energy, so this results in large fluctuations in the pressure. An alternative 
point of view is to realize that two phases are in equilibrium at the same 
temperature and pressure but with different energy densities, hlcncc at fixed 
P and T there can be large fluctuations in energy density. As in the case of 
the diatomic gas, whenever it is easy to transfer energy between collective 
and internal degrees of freedom, one should expect an enhancement of the 
bulk viscosity. The delay in attaining equilibrium due to energy exchange is 
ultimately responsible for entropy production. 

Generally speaking, one ought to expect a minimum in the shear viscosity 
and a maximum in the bulk viscosity near a phase transition. This is not 
a rigorous result but seems to be supported by various observations and 
calculations. The dimensionless ratio of shear or bulk viscosity to entropy 
(disorder) would seem to be a good way to characterize the intrinsic ability of 
a substance to relax towards equilibrium independent of the actual physical 
conditions (gradients of pressure, energy density, etc.). It would also seem 
to be a good way to compare very different substances. Unfortunately, the 
entropy density is not always measured or calculated, or is not made available 
to the reader, at the same time as the viscosities. This would facilitate more 
comparisons between atomic and molecular systems and strongly interacting 
matter at high temperature. 

3 Relativistic Dissipative Fluid Dynamics 

The energy-momentum tensor density for a perfect fluid (which does not 
imply that the matter is non-interacting) is T^*^ = —Pg^^ + wu'^u'^. Here 
w = P + e = Ts is the local enthalpy density and is the local flow velocity. 
Corrections to this expression are proportional to first derivatives of the local 
quantities whose coefficients are the shear viscosity 77, bulk viscosity and 
thermal conductivity x- Explicit expressions may be found in textbooks 
[27, 28] which are useful to summarize here. Dissipative contributions are 
added to the energy-momentum tensor and baryon current as follows. 
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= UBuf' + AJ^ (10) 

There arc two common definitions of the flow velocity in relativistic dissipa- 
tive fluid dynamics which are important to distinguish. 

In the Eckart approach is the velocity of baryon number flow. The 
dissipative terms must satisfy the conditions AJg = and u^u^AT'^'' = 0, 
the latter following from the requirement that T^^ be the energy density in 
the local (baryon) rest frame. The most general form of /ST'^^ is 

(11) 

Here 

H"'' = ui'u'' - g""" (12) 
is a projection tensor normal to u^^, 

A^ = d^-uydp (13) 

is a derivative normal to u^, and 

Qa = daT - Tu^dpUa (14) 

is the heat flow vector whose nonrelativistic limit is Q = — VT. The entropy 
current is 

= su>' + l;U^AT>''' . (15) 



T 



Its divergence is 



2T 

+ ^{V-nf + ^{VT + Tix)\ (16) 

All three dissipation coefficients must be non-negative to insure that entropy 
can never decrease. 

In the Landau-Lifshitz approach -u^ is the velocity of energy transport. 
The dissipative part of the energy-momentum tensor satisfles u^AT^^ = 0, 
and AJ^ is not constrained to be zero. In this case the most general form of 
the energy-momentum tensor is 

AT^"" = 77 {A^'u'' + A'^u^') + (§77 - C) H'^'dpuf . (17) 
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The baryon current is modified to 



A.g^x(^)^A.(^^). (18) 

Even though the entropy current in this approach is different, being 

= su^^ - ^AJ^ , (19) 

its divergence is the same. Physical, observable results cannot depend on 
how one defines the frames of reference. 

In nuclear collisions at RHIC the ratio of the net baryon number to the 
total number of particles is very small, typically less than 1%. With a greater 
energy available at the LHC it will be smaller yet. Typically hydrodynamic 
calculations set the net baryon number to zero as a good approximation. In 
that case one has no choice but to use the Landau- Lifshitz definition of fiow 
velocity. Then thermal conductivity is neither relevant nor defined. 

First order dissipative fluid dynamics applies when the viscosities are 
small, or when the gradients are small, or both. The dispersion relations for 
the transverse and longitudinal (pressure) parts of the momentum density 
are 

u + iDtk"^ = 
cu^ - v^^k'^ + iDi(jjk^ = (20) 

where — rj/w and Di — {^ri-\-C,)/w are diffusion constants with the dimen- 
sion of length and Vs is the speed of sound. Since w = Ts, the dimensionless 
ratio of shear or bulk viscosity to entropy density is a good way to character- 
ize the strongly interacting matter created at RHIC or LHC and to compare 
it to ordinary atomic and molecular systems. 



4 Shear Viscosity 

What are the theoretical predictions for the shear viscosity of hadrons and 
quark-gluon plasma? In the low energy chiral limit for pions the cross sec- 
tion is proportional to s//^, where s is the usual Mandelstam variable for 
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invariant mass- squared and f^^ is the pion decay constant. The thermally av- 
eraged cross section is (cr) oc T^/ f^, which leads to rj/s (x {f-^/T)"^. Explicit 
calculation gives [29] 

^ - —4 (21) 

Thus the ratio rj/s diverges as T ^ 0. At the other extreme lies quark-gluon 
plasma. The parton cross section behaves as cr oc g^/s. A first estimate yields 
r)/s cx 1/g'^. Asymptotic freedom at one loop order gives oc l/ln(T/AT) 
where 

At is proportional to the scale parameter Aqcd of QCD. Therefore rj/s is 
an increasing function of T in the quark-gluon phase. As a consequence, rj/s 
must have a minimum. Based on atomic and molecular data, this minimum 
should lie at the critical temperature if there is one, otherwise at or near the 
rapid crossover temperature. 

The most accurate and detailed calculation of the viscosity in the low 
temperature hadron phase was performed in [29]. (See also [30].) The two- 
body interactions used went beyond the chiral approximation, and included 
intermediate resonances such as the p-meson. To calculate the entropy den- 
sity a free gas consisting of pions, r) mesons, kaons, p mesons, cu mesons, and 
K* mesons was used. The results are displayed in figure 7, both two flavors 
(no kaons or K*) and three flavors (with kaons and K*). The qualitative 
behavior is the same as in eq. (21). 

The most accurate and detailed calculation of the shear viscosity in the 
high temperature quark-gluon phase was performed in [31]. They used per- 
turbative QCD to calculate the full leading-order expression, including sum- 
mation of the Coulomb logarithms. For three flavors of massless quarks the 
result is 

ri 5.12 , , 

- = (22) 

Use this together with the two-loop renormalization group expression for the 
running coupling 

1. 9 Inf^U 4 ln(21nf^)) (23) 
g^{T) 87r2 \AtJ \ \AtJJ ^ ' 

with At = 30 MeV, which approximately corresponds to using an energy 
scale of 2'kT and A^g^- = 200 MeV. The result is also plotted in figure 7 [19]. 
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These results imply a minimum in the neighborhood of the expected value 
of Tc ~ 190 MeV. Whether there is a discontinuity or a smooth crossover 
cannot be decided since both calculations are unreliable near T^. 

It is interesting to ask what happens in the large Nf. limit with g'^N^ 
held fixed [32]. In this limit, meson masses do not change very much but 
baryon masses scale proportional to N^,; therefore, baryons may be neglected 
in comparison to mesons due to the Boltzmann factor. Since the meson 
spectrum is essentially unchanged with increasing Nc, so is the Hagedorn 
temperature. The critical temperature to go from hadrons to quarks and 
gluons is very close to the Hagedorn temperature, so that Tc is not expected 
to change very much either. In the large Nc limit the meson-meson cross 
section scales as 1/iV^. According to our earlier discussion on the classical 
theory of gases, this implies that the ratio rj/s in the hadronic phase scales 
as A^"^. This general result is obeyed by (21) since it is known that scales 
as Nc- The large Nc limit of the viscosity in the quark and gluon phase may 
be inferred from the calculations of [31] to be 



where r = Nf/Nc- Thus the ratio rj/s has a finite large Nc limit in the 
quark and gluon phase. Therefore, we conclude that rj/s has a discontinuity 
proportional to N^ if Ac — ^ oo. This jump is in the opposite direction to 

that in figure 7. 

There are other means to compute the shear viscosity. Why not just 
compute it from numerical calculations with lattice gauge theory? The basic 
problem, which has been known for 25 years, is that one must determine the 
spectral density in frequency space from a Fourier transform of a function in 
(imaginary) time whose values are only known at a finite number of steps, 
typically A^^ = 4, 8, 16, 32, and then find the zero frequency limit. For exam- 
ple, quantitative results for the shear viscosity in lattice gauge theory have 
been reported by Nakamura and Sakai [33] for pure SU(3) without quarks. 
(This theory has a first order phase transition.) This bold effort obtained 
r]/s ^ 0.5 in the temperature range 1.6 < T/Tc < 2.2, albeit with uncer- 
tainties of order 100%. An improved calculation by Meyer [34, 35] obtained 
ri/s ^ 0.134 at T/Tc = 1.65 and r]/s ^ 0.102 at T/T^ = 1.24. Taken together 
these suggest a minimum at or near Tc. 




(24) 
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Another interesting approach has been the computation oi rj/s m M = 

4 supersymmetric SU{Nc) Yang-Mills theory (SYM). In that theory it is 
possible to do calculations in the large coupling limit [36] and in the weak 
coupling limit [37] , and without too much imagination it is possible to find a 
paramtrization that interpolates smoothly between the two limits. However, 
SYM has no renormalization group running coupling, no asymptotic freedom, 
and no thermodynamic phase transition. Since it has so many more degrees 
of freedom than QCD as possible scattering targets, its viscosity to entropy 
ratio is much smaller than QCD at high temperature when compared at the 
same value of the gauge coupling. However, when the theories are compared 
at the same value of the Debye screening mass they do agree reasonably well 
[37]. 

In another approach, Gelman, Shuryak and Zahed [38] have modeled the 
dynamics of long wavelength modes of QCD at temperatures from to 
1.5Tc as a classical, nonrelativistic gas of massive quasi-particles with color 
charges. They obtained a ratio of rj/s ^ 0.34 in this temperature range. 
One can also perform hadronic molecular dynamic simulations of hadronic 
collisions within a box to compute the shear viscosity from eq. (9). This has 
been done by Muronga [39] using the well-known UrQMD (Ultra-relativistic 
Quantum Molecular Dynamics) program which has been used so successfully 
for heavy ion collisions at energies below RHIC and as an afterburner in 
RHIC collisions. The practical problem with this approach is that reactions 
of the type 2 — > n particles are fairly well-known, and can be implemented 
in the computer simulation, but how to handle the inverse process n — > 2 
particles is not. Therefore all reactions involving n > 2 had to be turned off, 
otherwise detailed balance could not be achieved. The results were within 
a factor of 2 to 4 compared to those of Prakash et al. [29], depending on 
temperature. 

5 Bulk Viscosity 

What are the theoretical predictions for the bulk viscosity of hadrons and 
quark-gluon plasma? In the low energy chiral limit for pions the equation 
of state is e = 3P and = 1/3. Hence the bulk viscosity vanishes in the 
low temperature limit. However, there are deviations from this conformal 
behavior due to interactions among the pions which involve a characteristic 
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scale Ap ^ 275 MeV. Explicit calculation gives [40] 



c 



s 



87r2 
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(25) 



The ratio (/s goes to zero as T — > 0, in contrast io -q/s which diverges. 

At the other extreme lies quark-gluon plasma. At asymptotically high 
temperatures, much greater than quark masses and the QCD scale parame- 
ter A;^, the system is again essentially conformally invariant and the bulk 
viscosity should go to zero. The most accurate and detailed calculation of 
the bulk viscosity in the high temperature quark-gluon phase was performed 
in [41]. For three flavors of massless quarks the result is 



There are two significant differences compared to the corresponding expres- 
sion for the shear viscosity. First, the overall coefficient is four orders of mag- 
nitude smaller, and second, the factor g'^ appears in the numerator rather 
than in the denominator. The ratio C,/s goes to zero at high temperature 
because the strong coupling g^(T) goes to zero. Conversely it becomes large 
at low temperatures whose scale is set by A;^. 

The picture painted by massive versus massless pions is somewhat differ- 
ent. Figure 9 shows the results from [29] (see also [30]), which calculated the 
bulk viscosity from the elastic scattering of massive pions using experimen- 
tally inferred cross sections. As for the shear viscosity, the entropy density 
was calculated on the basis of a free gas of pions and 1], p and u mesons. 
Contrary to massless pions, the ratio C/s is a. decreasing function of T rather 
than an increasing function of T. From figure 9 one might conclude that 
C/s is a monotonically decreasing function of T even as one goes from the 
hadronic phase to the plasma phase. However, the figure is also consistent 
with a peak in the vicinity of 150 to 200 MeV. 

Indeed, other approaches strongly suggest a peak in (/s in the vicinity 
of the phase transition or rapid crossover. Using a combination of an exact 
sum rule derived from the low-energy theorems of broken scale invariance 
in QCD, lattice data on the equation of state which sees a large increase 
in e — 3P in the vicinity of Tc, and an ansatz for the correlation function 
of the energy-momentum tensor, Kharzeev and collaborators have estimated 



s 50001n(6.34/^) ■ 



(26) 
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the behavior of (/s above T,. [42, 43]. Their result is that (/s goes to zero 
rapidly at high temperature but becomes quite large as Tc is approached from 
above. For example, C/s < 0.01 for T/T^ > 1.4, (/s ^ 0.025 at T/T^ = 1.2, 
C/s ^ 0.07 at T/Tc = 1.1, and (/s ^ 0.3 at T/T^ = 1.01. The large T limit is 
numerically consistent with the perturbative QCD results displayed in figures 
8 and 9, and shows a very strong growth as T^. is approached. However, it 
should be mentioned that Moore and Saremi have challenged this analysis 



A direct calculation of C/s has been done by Meyer [45, 35] using lattice 
methods. Such a direct approach avoids the necessity to make an ansatz for 
the spectral density, but on the other hand it is difficult to access the low 
frequency limit of the spectral density. The results for pure gauge theory 
without quarks (this theory has a first order phase transition) are as follows: 
C/s < 0.015 at T/Te = 3.2, C/s ^ 0.008 at T/T^ = 1.65, and C/s f« 0.065 
at T/Tc = 1.24. The latter two points were obtained with a temporal lattice 
size of Nr = 12. For a smaller lattice of Nr = 8, the result was C/s ~ 0.73 at 
T/Tc = 1.02. Taking into account systematic and statistical uncertainties, 
both of these approaches agree with each other, from which we conclude 
that a sharp rise is predicted by lattice gauge theory as Tc is approached 
from above. 

Just as the ability to exchange energy between translational and internal 
degrees of freedom increases the bulk viscosity, so does the ability to store 
energy in, and extract energy from, a mean field. A simple formula was 
derived by Paech and Pratt [46] for the ubiquitous example of a cr field. 



Here the derivative of the pressure with respect to the mean field is carried 
out at fixed energy density. The damping coefficient is obtained from the 
dispersion relation of the field at zero momentum. 



Whenever the field a is rapidly varying or when rricr is small or both, such as 
near a phase transition, the bulk viscosity is predicted to be large. Within a 
simple mean-field model they found the behavior near the critical tempera- 



[44]. 




(27) 



(28) 
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ture of a first-order phase transition to be 



c 



—n 



(29) 



s 



with n = 1. In a more sophisticated approach the critical exponent is ex- 
pected to be somewhat different. Whether the ratio C/s diverges or not in 
full QCD is an interesting and open question [43]. 



The AdS/CFT (Anti-de Sitter/ Conformal Field Theory) correspondence, or 
more generally gauge/gravity duality, may offer new and profound insights 
into the non-perturbativc realm of QCD [47, 48, 49]. The bottom- up ap- 
proach, referred to as AdS/QCD, constructs a dual gravity theory with just 
a few parameters which are fit to some set of hadronic data. Then all other 
observables can be predicted. This approach has been of particular inter- 
est for calculating transport coefficients ever since the influential paper of 
Kovtun, Son and Starincts [13]. The basic ideas and results are summarized 
here; for a more extensive introduction to this subject, see [50]. 

Consider a background metric in 5-dimensional space of the generic form 



where the three functions goo, Qxx and Qrr depend on the fifth dimensional 
variable r only. The convention is that goo — —1 and g^x = 1 in flat 4- 
dimensional space. Our normal 4-dimensional world exists at r — > oo. There 
is a black brane (in four dimensions it would be a black hole) horizon located 
at r = To- As the horizon is approached it is assumed that the metric behaves 
as follows: grr — ^ 7r/ (^ — ^o) and goo — tq) while gxx approaches the 

finite value gxxifo)- The Hawking temperature is 



It is this black brane that provides for a finite temperature system. The 
matter and radiation fields which generate this metric are oftentimes not 
known. The procedure then is to compute the correlation functions of the 



6 Gauge/Gravity Correspondence 



ds^ = goo{r)dt'^ + gxx{r)\dx.f + grr{r)dr'^ 



(30) 




(31) 
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energy-momentum tensor to infer the long wavelength dispersion relations 
from (i) the poles of these functions or (ii) the Kubo formulas, (iii) to use 
the membrane paradigm where one solves the field equations with boundary 
conditions applied on a stretched horizon of the black brane, or (iv) to solve 
for gauge invariant perturbations with approriate conditions. The results 
from all of these methods coincide, as demonstrated for cases (i), (ii) and 
(iv) in [51]. See [52] for case (iii). 

A quite general formula for the ratio of shear viscosity to entropy density 
was first derived in [53] and reproduced by different methods in [54]. 



!L = T^J^E^^ Tdr^p^kl^. (32) 
^ y-9oo{ro)grr{ro) •^'■0 \J -g{r)gxx{r) 

Hence, knowing the background metric one may compute rj/s just by eval- 
uating this integral. A frequently used metric is one that is dual to jV" = 4 
SU(A^) supersymmetric Yang Mills theory at finite temperature and in the 
limit N ^ oo and g'^N — > oo. 



l2 

-f{r)dt^ + \d^\^ +^rjrTdr^. (33) 
J r^/(r) 



Here /(r) = 1 — (ro/r)^. Using this metric yields rj/s = l/47r, which is 
the famous result. In fact, there is a theorem which states that whenever 
ii!o = R%, where R^i, is the Ricci tensor, then rj/s = l/in [55]. This equality 
of components of the Ricci tensor is true whenever the background metric is 
generated by scalar fields [56]. 

The gravity dual to QCD is not known; in fact, it is not even known if 
one exists. AdS/CFT can only be an approximate representation of QCD 
at very high energies or temperatures. The reason is that a conformal field 
theory has no intrinsic energy scale, although it may have dimensionless 
parameters. This is approximately representative of QCD at high energies 
or temperatures since the QCD gauge coupling varies only logarithmically. 
There are several phenomenological models for QCD that incorporate con- 
finement. The simplest of these is the hard wall model [57] which uses the 
metric of Eq. (33) but with r-space cutoff at some value r„^i^. This leads 
to radial excitations of the vector and axial-vector meson spectra with the 
mass being linear in the radial quantum number n. Introduction of the new 
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scale rmin implies that the transport coefficients might vary nontrivially with 
temperature T. 

The soft wall model was developed to improve upon the hard wall model. 
In particular, it leads to linear Regge trajectories wherein the radial excita- 
tions of the vector and axial-vector meson spectra have mass-squared being 
linear in the radial quantum number n, in substantial agreement with data. 
It can be obtained by adding a dilaton field to the usual AdSs metric to 
break the scale invariance. The dilaton profile which leads to Regge behav- 
ior is 0(r) = cL^/r^, where c is a constant which can be determined by fitting 
the meson spectrum. If one computes the ratio ri/s using the soft wall model 
of [58] and Eq. (32), the result is 1/4tt because the metric is exactly AdSr,. 
There is an alternative formulation of the soft wall. Instead of adding a 
nontrivial dilaton, one keeps the dilaton constant while deforming the metric 
away from AdSs. This version of the soft wall model has been studied in 
[59]. The deformed metric is 



(34) 



It should be stressed that this is a different implementation of the soft wall 
model than that of [58]. Both pictures lead to linear Regge trajectories, but 
other physical quantities may differ. Using this metric in Eq. (32) yields 



7] 



s 27:xq 

where 



l + l(e--l) 



(35) 



xo = ^ . (36) 

The Hawking-Pagc analysis of the phase transtition in [60] gives Tc = l/(2y^) 
nip/ A ^192 MeV. Note that the conformal limit of l/47r is approached from 
below as T increases. One may be surprised that the shear viscosity is not 
greater than or equal to l/47r. The reason is that the theorem of [55] is not 
applicable because EP^^ R^.. Instead, in this version of the soft wall model, 
one finds by explicit computation that 

i?°-i?: = -60'(r)e2^W^ (37) 
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which is nonzero as long as c6(r) is not constant. This may be a sign that a 
metric such as (34) cannot be generated by conventional supergravity matter 
fields. It remains an open question. 

As discussed earlier, the bulk viscosity vanishes for conformally invariant 
theories, in particular AdS/CFT. So far, there is no known formula for the 
bulk viscosity analogous to Eq. (32). It seems necessary to know what fields 
generate the background metric. The simplest case is when the metric is 
created by a single scalar field with action 

s = / d'xV^ [r - \dyd^(t> - y (0)] (38) 

where is the gravitational constant in five dimensions, R is the curvature, 
and V(0) is a potential. Examples of such potentials are given by Chamblin 
and Reall [61] 

V{4>) = V^e'l' (39) 

and by Gubser et al. [62] 

y ((/)) = - ^ cosh (70) + 60' (40) 

where Vq, 7, and h are constants. These and related potentials have parame- 
ters that can be varied to approximate the QCD equation of state as obtained 
from lattice QCD calculations. For example, for the Chambfin-Reall poten- 
tial the speed of sound is 

In the absence of a top-down approach it is very difficult to know which 
potential provides the best representation of QCD. 

Buchel [63] has conjectured that there is a minimum value of the bulk 
viscosity relative to the shear viscosity on the basis of a variety of holograph- 
ically dual computations. This inequality is 

^ > 2 (i - vi) . (42) 

Rather than become engrossed in detailed models, let us take a more general 
point of view to see what can be learned. 
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Motivated by the above studies, Springer [56] considered the following 
constraint on the metric 

goo{r) = -^gUr) + a,gll-\r) . (43) 
The requirement that R^q = R% leads to 

"'"'If'''' {i^-lMrnJ. (44) 

The parameter oi is not independent but is fixed in terms of the others by 
the requirement that Qrr — 7r/('" — '"o) and 5^00 — 7o('" — '"o) while Qxx 
approaches the finite value QxxiT^o)- The Hawking temperature is 

T-^glx^ro). (45) 

The speed of sound and the viscosities can be found by calculating the dis- 
persion relation. 

2 00 — 02 

77 1 



s An 
s Stt 



(46) 



These results respect the conjectured minimum value of C/^ if + 02 > 3; 
the equality is attained in certain models. By specifying the metric one 
may deduce the speed of sound (hence the equation of state) and the bulk 
viscosity. The advantage of this approach is that it allows one to correlate 
the equation of state and the dissipative coefficients in a consistent way. The 
disadvantage is that until we have an approximate or exact gravity dual to 
QCD, this approach will remain phenomenological. 



7 Observable Consequences 

It ought to be possible to extract numerical values of the viscosities in heavy 
ion collisions via scaling violations to perfect fluid flow predictions [64] . The 
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program is to solve relativistic viscous fluid equations, with appropriate ini- 
tial conditions and with a hadron cascade afterburner [65], over a range of 
beam energies and nuclei and extract r}{T)/s{T) and ({T)/s{T) from com- 
parison with data. This program is analogous to what was accomplished at 
lower energies of 30 to 1000 MeV per nucleon beam energies in the lab frame. 
At those energies, it was possible to infer the compressibility of nuclear mat- 
ter and the momentum-dependence of the nuclear optical potential via the 
transverse momentum distribution relative to the reaction plane [66] and via 
the l)alaucc between attractive and repulsive scattering [67]. 

To make quantitative comparisons with data it is necessary to solve rela- 
tivistic viscous fluid equations numerically in three space dimensions without 
assuming any particular symmetry. On general grounds one should expect 
viscous effects to smooth out gradients in temperature and flow velocity and 
to slow down the expansion of the system. To develop some intuition as to 
what happens it is worthwhile to solve those equations in one space and one 
time direction under the assumption of boost invariance [68]. This is the 
fluid dynamical model proposed by Bjorken [69]. 

For the sake of simplicity, consider a bag model type equation of state at 
temperatures above the critical one. The pressure is 



and the energy and entropy densities follow directly. Neglect the shear vis- 
cosity in order to focus attention on the bulk viscosity. The fluid dynamic 
equation to solve is 



where t denotes the local time (we do not use the symbol r as is customarily 
done so as not to confuse it with the relaxation time discussed earher) . First 
make the assumption that C is a monotonically decreasing function of T to 
the inverse power of n > 0. 



Here the subscript i indicates the initial value at the beginning of the fluid 




(47) 




(49) 
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expansion. The solution is easily found. 



n + A Q 1 
n + lsi UTi 



The second term in the square bracket should be small comapared to one in 
order that the gradients are not so large initially as to negate the application 
of viscous fluid dynamics. The last term on the right hand side shows that 
the temperature decreases with time more slowly in the presence of the bulk 
viscosity. In fact that term dominates at large t. Eventually the temperature 
will fall to Tc and one would need to switch to the hadronic equation of state 
and associated viscosity. 

Next make the assumption that C has a power singularity at Tc- 

cm-c.(^)" (51) 

The equation of motion cannot be solved in terms of elementary functions. 
However, the asymptotic behavior can be determined to be 




r(i)^Te+(Ti-Te)P — ast^oo. (52) 



It takes an infinite amount of time to reach the critical temperature. This is 
an example of critical slowing down. Of course this would not happen in the 
real world, because eq. (48) neglects transverse expansion and the surface 
emission of particles near the cooler periphery of the matter. 

At RHIC and LHC some of the specific proposals to extract or infer 
the viscosity to entropy ratio from data include: elliptic fiow [70] , Hanbury- 
Brown and Twiss (HBT) interfcromctry [70], single particle momentum spec- 
tra [70, 71], and momentum fluctuations [72]. Other possibilities include jet 
quenching and photon and dilepton spectra. Some of the complications in- 
clude the possibility that gradients are so large that the second-order dissi- 
pative equations of Israel and Stewart are necessary [73] and that turbulence 
in the plasma may lead to an anomalous viscosity [74]. Clearly this is an 
interesting and challenging goal but worth the effort. 
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8 Conclusion 



There are strong arguments and calculations which suggest that hadron/quark- 

gluon matter should have a minimum in the ratio of shear viscosity to entropy, 
and a maximum in the ratio of bulk viscosity to entropy, at or near the crit- 
ical or crossover point in the phase diagram. Suffiently detailed calculations 
and measurements ought to allow us to infer these quantities from experi- 
ments at RHIC and LHC. These are interesting dimensionless measures of 
dissipation relative to disorder which can be compared among a wide variety 
of substances. 

My conclusion is that RHIC is a thermometer since it measures hadron 
ratios and photon and dilepton spectra, it is a barometer since it measure 
radial and elliptic flow, and it may be a viscometer since it could measure 
deviations from ideal (nonviscous) fluid flow. Soon the same will happen 
with the LHC. There is much work still to be done by both theorists and 
experiment alist s ! 
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Figure 1: The ratio 77/s as a function of T for helium with s normahzed 
such that s(T = 0) = 0. The curves correspond to fixed pressures, one of 
them being the critical pressure, and the others being greater (1 MPa) and 
the other smaller (0.1 MPa). Below the critical pressure there is a jump in 
the ratio, and above the critical pressure there is only a broad minimum. 
They were constructed using data from NIST and CODATA. Prom Csernai, 
Kapusta and McLerran. 
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Figure 2: The ratio 77/5 as a function of T for nitrogen with s normahzed 
such that s{T = 0) = 0. The curves correspond to fixed pressures, one of 
them being the critical pressure, and the others being greater (10 MPa) and 
the other smaller (0.1 MPa). Below the critical pressure there is a jump in 
the ratio, and above the critical pressure there is only a broad minimum. 
They were constructed using data from NIST and CODATA. The curves are 
plotted on logarithmic scale to make the behavior around the critical point 
more visible. Prom Csernai, Kapusta and McLerran. 
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Figure 3: The ratio r]/s as a function of T for water with s normahzed 
such that s(T = 0) = 0. The curves correspond to fixed pressures, one of 
them being the critical pressure, and the others being greater (100 MPa) and 
the other smaller (10 MPa). Below the critical pressure there is a jump in 
the ratio, and above the critical pressure there is only a broad minimum. 
They were constructed using data from NIST and CODATA. Prom Csernai, 
Kapusta and McLerran. 
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Figure 4: Scaled ratios along the critical isobar for helium, nitrogen and water 
obtained from the previous three figures. There is no obvious universality. 
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Figure 5: The shear viscosity versus the Coulomb couphng parameter for a 
two-dimensional Yukawa system in the liquid state. The minimum in the 
viscosity occurs when it has approximately equal contributions from kinetic 
and potential components. Prom Liu and Goree. 
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Figure 6: The bulk viscosity versus density at fixed temperature for a clas- 
sical liquid-gas phase transition using a Lennard- Jones potential. The bulk 
viscosity has a peak near the critical point corresponding to the temperature 
T* — 1.35. Prom Meier, Laesecke and Kabelac. 
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Figure 7: The ratio r}/s for the low temperature hadronic phase and for the 
high temperature quark-gluon phase. Neither calculation is very reliable in 
the vicinity of the critical or rapid crossover temperature. From Csernai, 
Kapusta and McLerran. 
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Figure 8: The ratio C/s for the low temperature hadronic phase represented 
by chiral pions and for the high temperature quark-gluon phase. Neither 
calculation is very reliable in the vicinity of the critical or rapid crossover 
temperature. 
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Figure 9: The ratio C/s for the low temperature hadronic phase represented 
by massive pions and for the high temperature quark-gluon phase. Neither 
calculation is very reliable in the vicinity of the critical or rapid crossover 
temperature. 
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